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Abstract 

The scaling behavior of the mutifractahty for the sea-bottom topography 
of the South Sea in Korea is numericahy investigated. In particular, we focus 
on the behavior of the gth-moment depth-depth correlation function of the 
sea-bottom topography and its multifractal spectrum. Through the multi- 
fractal analysis, the fractal dimension and the scaling exponents are obtained 
numerically, and the relation between the Hurst exponent and the fractal di- 
mension is also derived. 
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For more than two decades, there has been considerable interest in the in- 
vestigation of the scaUng behavior on fractal models. The pioneering work[l] 
on fractals has introduced the concept of fractals and shown some relation 
between the self-similar fractals and the self-affine fractals. The self-afHne 
fractals [2 — 5] constitute random and complicated structure that has been 
applied to broader range of problems such as Eden and ballistic deposition 
model[6 — 9], mountain heights, clouds, coast lines and cracks. Among other 
examples of many fractal models we also can mention the self-avoiding random 
walk, percolation clusters, diffusion-limited aggregations, random resistor net- 
works, polymer bonds, turbulences, chaotic motions [2,4, 5, 10 — 13], etc. 

On the other hand, Matsushita and Ouchi[14] have recently shown that 
the self-afHne of mountain topography in Japan is mainly obtained by using 
the numerical method, from which they have discussed the relation Df = 
2 — H = 1.37 between the sclf-affine exponent H and the fractal dimension 
Df for the standard deviation of mountain heights. More recently, Katsuragi 
and Honjo[15] have dealt with self-affine fractal profiles on the entropy spec- 
tra method and for the whole range of a mountain they have obtained the 
generalized multi-affine profiles and the entropy spectra. 

It has recently discussed that the multifractals[2, 5] have essentially the 
distributions of singularities of the scaling exponents on interwoven sets of 
various fractal dimensions. The multifractal formalism for the box-counting 
method is represented in terms of their generalized dimensions and spectra. 
The multifractality of self-affine fractals has been discussed by Barabdsi et 
aZ[16, 17]. They have investigated the multi-affine function and the multifrac- 
tal spectra, and calculated the generalized Hurst exponent from the g-th order 
height-height correlation function. Very recently, the generalized dimension 
and the scaling exponent which are interesting in multifractals for mountain 
heights in Korea have been calculated by Kim and Kong[18]. We think that 
it is interesting to be extended to the multifractals, although we do not have 
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enough sets of data for the seabottom depth. In particular, the multifractal 
theory is a good basis for analyzing the statistical data on the seabottom 
depth, and really provides the probabilistic evidence for the relation between 
Hurst exponent and the fractal dimension. 

The main purpose of this paper is to investigate numerically the behav- 
ior of the depth-depth correlation function and its multifractal spectrum in 
the seabottom depth from the seabottom topography between the longitude 
129°21' - 129°54'£; and the latitude 34°52' - 35°19'iV located near to the 
east of Pusan. The primary reason we select this region is that this region 
has more data for sea-bottom depth than those of any other region around 
Korean peninsula. We consider well known relations of the depth-depth cor- 
relation function and the multifractals. The q'th-moment depth-depth correla- 
tion function is also numerically discussed from the multi-scaling properties of 
seabottom depth function. In addition, from the extrapolated sets of data for 
the seabottom depths, we briefly explain an efficient and convenient method 
relevant to numerical investigation on the scaling behavior of the generalized 
dimension and the spectra. 

First of all, we are interested in the scaling quantities obtained from the 
extrapolated data of the seabottom depth on a two dimensional square lat- 
tice between 129°21' - 129°5A'E and 34°52' - 35°19'N, as shown in Fig. 1. 
Moreover, the current data are taken for the values of the seabottom depths 
projected on 100 x 100 lattice points in square area of 50 x 50 km^, where the 
distance between one lattice point and its neighboring lattice point is 500 m. 
Let's consider the gth-moment depth-depth correlation function Cq{r) defined 

by 

1 ^ 

C,{r) = -J2\D{ri + r)-D{riW, (1) 

i=l 

where D{ri) is the ith depth in the r-direction and r denotes two perpendic- 
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ular longitudinal or latitudinal direction on a two dimensional square lattice. 
This correlation function can be described as a non-trivial multi-scaling be- 
havior; 

C,(r)~r'^^^ (2) 



where Hq is the generalized Hurst exponent in the limit of r — > 0. 

Next, we discuss the generalized dimension Dq, the scaling exponents aq, 
and fq on the multifractal structure. In general, it has been known that the 
generalized dimension is represented as the fractal distribution having the sin- 
gular vahics infinitely. Accordingly, the generalized dimension and the scaling 
exponents are known to be given by following equations[4, 5] 



1 In^^ (3) 
^ e-^oq - 1 Ine 



1 Y^iPi^lnpi 



and 



€- 



where pi is the probability of the seabottom depth existing on the i-th box with 
the square area of e x e and the scaling quantity is the number of the box 
having the probability pj. By introducing the above expressions, the spectra 
fq and aq are simply calculated from eqs.(4) and (5) with fq = qaq — {q — l)Dq 
and aq = ■^[(q—l)Dq], where fq and Oq can be obtained by Legendre transfor- 
mation. In our scheme, we will make use of eqs.(3)-(5) to find out the fractal 
dimension and other scaling exponents, and these mathematical techniques 
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lead us to more general results. 

We present more detailed numerical data of the depth-depth correlation 
function of the seabottom depth and its multifractal spectrum. As mentioned 
above, we assume that the seabottom depths divided by the intervals of 500 
m are located on each region of two dimensional lattice. As listed in Table 

1, wc take into account the data of three regions ^4 - C in seabottom depths 
between 110 m and 125 m. As we only restrict ourselves to the longitudinal 
direction from our data, the generalized Hurst exponent from the gth-moment 
depth-depth correlation function is calculated numerically. Hence, from eq.(2) 
the generalized Hurst exponents H2, H^, and H4 on the longitudinal direction 
are found to be 0.599, 0.587, and 0.586, respectively. 

In our multifractal structure, the box-counting method is used for two 
square areas of 2.5 x 2.5 km^ and 5x5 km^. Through the multifractal anal- 
ysis, the generalized dimension and the scaling exponents are also obtained 
from the data with the number of seabottom depths in Table 1. These values 
by using the theoretical expressions of eqs.(3) — (5) are calculated numerically 
in our three regions. In Fig. 2 we show the scaling exponents and fg on 
the multifractal structure, and it is found numerically that the scaling expo- 
nents are estimated as a+oo = 1.097696 and CK—oo — 2.172183 only in the case 
of the region A. Particularly, as shown schematically in Fig. 3 and Table 

2, it is obtained that the maximum values of the generalized dimension, i.e., 
in our three regions A, B, and C the fractal dimension Df = Dq are respec- 
tively calculated as 1.312476, 1.366726, and 1.372243 for the case of the square 
area 2.5 x 2.5km?. Hence, we obtain that Hurst exponent takes approximately 
0.65 in our multifractal structure, while it takes the value H=0.63 in Japanese 
mountain topography [14]. 

In conclusion, we have remarkably studied on the gth-moment depth-depth 
correlation function and its multifractal spectrum in the seabottom topogra- 
phy of the South Sea in Korea, as shown in Fig.l. Specifically, the multifractals 
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have been investigated numerically by employing the box-counting method 
on a two dimensional square lattice. Based on these results, we expect that 
further analytical and numerical progress for multifractals may be achieved 
from more data for the ocean floor on diverse region of seabottom topogra- 
phy. Futhermore, it will be useful to apply and to reinvestigate the above 
formalism in many scientific ficlds[19, 20] such as quantum disorder systems, 
irreversible growth structures, and fractured surfaces of minerals. 

This work was supported in part by the academic research fund of Puky- 
ong National University in Korea. 
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Figure Captions 



Fig. 1 : The seabottom depth in seabottom topography taken from 3^2^ 
scale map between the longitude 129"21' - 129''54'^ and the latitude 34"52' 
- 35°19'A^, where the interval of seabottom depth between the solid lines is 
found to be 10m. 

Fig. 2 : fq versus aq by using the box - counting method for the case of 
the square area 2.5 x 2.5km^. The values of the seabottom depths in three 
regions A, B, and C are respectively given by the thin solid, dashed, and thick 
solid lines. 

Fig. 3 : Plots of aq and fq as a function of q on the seabottom depths for the 
case of the square area 2.5 x 2.bkm^, where the vertical bars are error bars 
averaged over the square areas of three regions A, B, and C. 



Table Captions 

Table 1 : Number of the seabottom depths in our three regions A - C between 
110m and 125m on the seabottom topography. 

Table 2 : Summary of values of Dq, aq, and fq calculated from the data 
of three regions A, B, C for the case of the square area 2.5 x 2.5A;m^. 
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This figure "figurel.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0006107v2 



This figure "figure2.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0006107v2 



This figure "tablel-2.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0006107v2 



This figure "figure3.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0006107v2 



